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Abstract 



> 

^o . 

(^ , We construct a Madelung fluid model with specific time variable 

VO I parameters as dissipative quantum fluid and linearize it in terms of 

lO ' Schrodinger equation with time dependent parameters. It allows us 

^^ , to find exact solutions of the nonlinear Madelung system in terms of 

solutions of the Schrodinger equation and the corresponding classi- 
cal linear ODE with variable frequency and damping. For the com- 
plex velocity field the Madelung system takes the form of a nonlinear 
r> I complex Schrodinger-Burgers equation, for which we obtain exact so- 

C^ ' lutions using complex Cole-Hopf transformation. In particular, we 

discuss and give exact results for nonlinear Madelung systems related 
with Caldirola-Kanai type dissipative harmonic oscillator. 

1 Introduction 

In the recent years the Madelung fluid description of quantum mechanics has 
been applied to some flelds where the quantum formalism is a useful tool for 
describing the evolution of classical (quantum-like) systems and studying the 



dispersionless or seiniclassical limit of nonlinear partial differential equations 
of Schrodinger type, pQ. The Madelung fluid representation, proposed first 
by [2], being a complex quantity, represents a solution of the Schrodinger 
equation, in terms of modulus and phase. Substituted to the Schrodinger 
equation it allows to obtain a pair of nonlinear hydrodynamic type equa- 
tions. Thus, the Madelung fluid equations are nonlinear system of PDEs, 
while the Schrodinger equation is the linear one. Then, the Madelung trans- 
form is a complex linearization transform, similar to the Cole-Hopf trans- 
formation, linearizing the nonlinear Burgers equation in terms of the linear 
heat equation, see |1] and [5]. Nonlinear models admitting such type of direct 
linearization are called by F. Calogero as C-integrable models. 

In this work, we construct a Madelung fluid model with time variable pa- 
rameters as dissipative quantum fluid and linearize it in terms of Schrodinger 
equation with time dependent parameters. It allows us to find exact solution 
of the nonlinear Madelung system in terms of solutions to the Schrodinger 
equation and the corresponding classical linear ODE with variable frequency 
and damping. Moreover, the Madelung system written for the complex veloc- 
ity field takes the form of a nonlinear complex Schrodinger-Burgers equation, 
which exact solutions we obtain using complex Cole-Hopf transformation. As 
known, in the usual Cole-Hopf transformation zeros of the linear heat equa- 
tion lead to poles in the corresponding Burgers equation. Similarly, in our 
case, by the complex Cole-Hopf transformation zeros of the Schrodinger equa- 
tion transform to pole singularities in the complex Schrodinger-Burgers equa- 
tion. Thus, using exact solutions of the linear problem, one can find also the 
dynamics of the poles in the corresponding nonlinear problem. As an exactly 
solvable model, we describe a dissipative nonlinear complex Schrodinger- 
Burgers equation of Caldirola-Kanai type, [TT], [12]. Exact solutions of the 
nonlinear models are found and the motion of zeros and poles is discussed 
explicitly. Some illustrative plots are constructed. 



2 The Schrodinger Equation and its Madelung 
Representation 

2.1 Solution of the Schrodinger Equation 

Consider the one-dimensional Schrodinger equation for harmonic oscillator 
with time-dependent parameters 



'dt ~ ~2fi{t) (9g2 



^^^ = -77^^ + ^^^^^^^^^^'*' (1) 

and initial condition 

^(g, to) = ^(g), -oo < g < oo. (2) 

Using the Evolution operator method, [8j, it was proved that, see [7], if x{t) 
is the solution of the classical equation of motion 

x + ^x + uj\t)x = 0, x{to)=xo^O, x{to) = 0, (3) 

fi{t) 

then the solution of the IVP dl])-© is found as ^(g,t) = f7(t, to)V^(g), where 
the evolution operator is 

U = exp (1/g^) exp (Mg| + 1)) exp (-1,|,) 

and the auxiliary functions are 



fi x(ty =■>' '"'J i^ioxm """ • " \x{t)\ 

In particular, if the initial function is the normalized eigenstate corresponding 
to eignenvalue Ek = h ^oik + 1/2) of the Hamiltonian for the standard 
harmonic oscillator, that is 

Mq) = Nke-^''"Hk{^oq), fc = 0,1,2,..., (4) 

then, the time-evolved state for the Schrodinger equation (^ is 

^fc(g,t) = ?7(t,to)</'fc(g) 



where 



Nk\jR{t) X exp ( M ^ + 7j ) arctan(r^o5'(^)) 
X exp (-^R^{t)q^\ X Hk (Jn'oR{t)q 



m 



Xn 



^x^{t) + {nox{t)g{t))\ 
The corresponding probabihty density is then 



(5) 



(6) 



Pkiq,t) 



2'=A;!0F 



X jQQR{t) X exp 



'noR{t)q 



^Hl 



iVLQR{t)q 



2.2 Madelung representation 

As known, Madelung representation of the complex-valued wave function 

^{.q.t) = Vpexp Qs) = exp (^-\np+^S^ , (7) 

where p = p{q,t) is the probability density and S = S{q,t) is the action, 
both being real-valued functions, decomposes the Schrodinger equation ([1]) 
into a system of nonlinear coupled partial differential equations. 



dS 1 

+ 



dt 2p{t) \dqj 



dSy ^ p{t)u\t) ^. 



h' 



2p{t) 



J_d\/p 
y/P dq"^ . 



dp d 
dt dq 



P 



1 dS 
pit) dq^ 



0. 



(8) 



The first equation may be viewed as a generalization of the usual Hamilton- 
Jacobi equation. The term with explicit Ti dependence is the quantum poten- 
tial encoding the quantum aspects of the theory. When ^ — )■ 0, the equation 
becomes Hamilton- Jacobi equation for a non-relativistic particle with time 
dependent mass. The second equation is a continuity equation expressing 



the conservation of probability density. Using the relation ([7]), one can see 
that the system ([H]), with general initial conditions 

S{q,to) = S{q), p{q,to) = p{q), (9) 

S{q),p{q) being real- valued functions, has formal solution 

S{q,t) = -zh\n (^1^) , PiQ,t) = m<l,t)\^ (10) 

where "^{qyt) is a solution of the Schrodinger equation ([T]) with initial con- 
dition 



^(g,io) = VP(g)exp(-5(g)). (11) 

We remark that, since ^I^(g,t) is complex-valued, in general S{q,t) is multi- 
valued, i.e. S{q, t) = -in ln(^/|^|) -f- 27m^, n = 0, ±1, ±2, ..., but fixing the 
initial condition S^q.to) = S{q) leads to a single-valued solution of the IVP. 

2.3 Madelung Hydrodynamic Equations 

1 dS 
Introducing classical velocity, v{q,t) = — rr^— , the system ([8]) transforms 

to Madelung fluid equations 



dv iiit) dv Id 



-n' ( 1 d'^p\ ^ p{t)u\t) ^. 



2p{t) Vv/p dq^ 



(12) 



dt p,{t) dq p{t) dq 

9p d , , ^ 

These equations are similar to the classical hydrodynamic equations where 
p{q,t) is the density and v{q,t) is the velocity field of the one-dimensional 
fluid. The system of fluid equations f lT^ with general initial conditions 

v{q,to) = v{q), p{q,to) = p{q), 

v{q),p{q) being real- valued functions, has formal solution 

where \l/(g, t) is solution of the Schrodinger equation ([T]) subject to the initial 
condition 



3 Potential Schrodinger-Burgers Equation 

Writing the wave function in the form 

^(g,t)=expQ/x(t)F(g,t) 

where F{q,t) is a complex potential, the IVP for the Schrodinger equation 
(II]) transforms to the following IVP for the nonlinear potential Schrodinger- 
Burgers equation 



'dt^'iKf) 2^W "^^"^ ~2^{t)dq^ ' (14) 

F(g,to) = F(g), 

Therefore, the formal solution of this problem is 

i^(9,t) = -^(ln*(g,t)), (15) 

where ^^(g, t) is solution of the Schrodinger equation ([T]), with initial condition 

^(g,to)=expQ/i(to)F(g) 

Note again, that fixing the initial condition we obtain a single- valued solution 
F{q,t). Now, using the Madelung representation (JTj) and relation (IT5l) . one 
can write 

ffe,) = F, + .F. = -ljS-^l„p. (16) 

where Fi = Fi(g, t) represents the velocity potential, and F2 = F2(g, t) the 
stream function of the fluid, ( Fi,F2 being real- valued) . Accordingly, the 
real and imaginary parts of the potential Schrodinger-Burgers equation (1141) 
become 



dt n 2 \ dq dq J 2 2/i dq'^ /-, y\ 

Of 2 Ac 9Fi9F2 h d^Fi ^ ' 



+ -F: 



dt /i dq dq 2fi dq^ 



Using the relations fll6p and fllOl) , one can see that the nonhnear system (II 7p 
with general initial conditions 

Fi(g,to) = Fi(g), F2(g,to) = F2(g) 

and Fi(g), F2(g) real-valued functions, has solution of the form 

where "^{q, t) is solution of the Schrodinger equation ([1]) with initial condition 
^(g,to) = exp f-/i(to)Fi(g)j x exp f --/i(to)F2(g) j . 

4 Schrodinger-Burgers Equation 

Representation of the wave function in the form 

^{q, t) = exp Q/x(t) J" Vi^, t)de) , (18) 

where V{q, t) is a complex velocity, transforms the IVP for the Schrodinger 
equation ([T]) to the following IVP for a nonlinear Schrodinger-Burgers equa- 
tion with time dependent coefficients 

dV fiit)^^ ^BV 2. N i^ d'^y 

■W + W*'^lk^^^'"^MftW (19) 

V(q,to)^V{q), 
Solution of this IVP is found by the complex Cole-Hopf transformation 

where '^{q,t) is solution of the the Schrodinger equation ([1]), corresponding 
to initial condition 

^(g, to) = Hq) = exp (^/i(to) /' ^(0^?^) • 



Using the Madelung representation ([7]) and the complex Cole-Hopf trans- 
formation fl2Up one can write the complex velocity function in the form 

Tw N . I dS ih d ,^ . , , 

VM^. + ^u^—--—-Hnp), (21) 

where v = v{q,t), u = u{q,t) are real- valued, v represents the classical 
velocity, and u the quantum velocity. This splits the Schrodinger-Burgers 
equation into real and imaginary parts, respectively. 



dv flit) dv 2/ \ ~^ '9^'" "^M 

du iiit) dv du h d'^v 



(22) 



dt fi{t) dq dq 2/i(t)ag2- 

The first equation is a hydrodynamic equation for the classical velocity, and 
the second one is the transport equation for the quantum velocity. Using re- 
lations fl2T|) and f lTOj) . we find that the system of nonlinear coupled equations 
( !22|) with general initial conditions 

v{q,to) = v{q), u{q,to) = u{q) 

has formal solution 

In kr^ , u = — — — in(l^l). 



where ^ = \E'(g, t) is a solution of the Schrodinger equation ([T]) with general 
initial condition 



ifx{to) /"^-/^x ,^x , fJ-it 



0) 



q 



5 Exactly Solvable Nonlinear Models 

The Caldirola-Kanai model, |Tl],[T2], which is a one dimensional system with 
an exponentially increasing mass, is the best known model of harmonic os- 
cillator with time-dependent parameters. Here, using the general discussion 
in the previous parts, we obtain exact solutions of the nonlinear problems 
related with the Caldirola-Kanai oscillator 

'^dt~ 2" dq' +2^°" ^^' (23) 

[ *(Q,0) = V(g), qeR, 

8 



where /i(t) = e'^ is the integrating factor, r(t) = 7 is the damping term, 
7 > 0, and (jj'^it) = Wq is a constant frequency. As known, solutions of 
the Cadirola-Kanai oscillator can be found in terms of the solution to the 
corresponding classical equation of motion 



X + 7i + WqX = 0, a;(0) = xq 7^ 0, i;(0) = 0. 



(24) 



Clearly, according to the sign of Vl? = oJq ~ (7^/4) there are three different 
type of behavior- critical damping, under damping and over damping. In 
this article, we discuss the critical damping case, i.e. Vl? = 0. If Q"^ = 
Uq — (7^/4) = 0, then the classical equation (^^ has solution 





x{t)-- 


= xoe ^(l + |t). 


and it follows that 






9it) = 


-nt 

i + it' 


/ pit 


^^' [ii + W + <t' 



1/2 



Then, using ([5]), exact solutions of the Schrodinger equation fl23|) with initial 
conditions \l/(g, 0) = ^Pkiq), are 



*fc(g,t) 



=7* 



1/4 



N. 



;i + ¥)' 






exp ii{k + -) arctan(— ^; 



M 







1- 



1 



(1 + |t)2 + wit 



2j-2 



o-yt 



'"'"'''^y 2h[il + Uy + w'ot\ 



l']H.{,^ 



=7* 



h \{i + w + wit 



2+2 



1/2 



a. Madelung representation of Caldirola-Kanai Oscillator. Madelung 
representation of the wave function decomposes the Schrodinger equation 
f l23|l into a system of nonlinear coupled partial differential equations. 



ds , 1 _,, /asV , 1^1 ,, 2 V 



dp 



7* 



d 



dt dq 



dS 



2 
0. 



1 d'^p 



This system of equations with specific initial conditions 



S{g, 0) = 0, ft(,, 0) = iV= exp ( -"^f] Hi ("#, 



has exact solutions of the form 



Skiq,t) 



UJn 



te<^ 



2 Vl + i^ 



1 



1 



[l + ltf + wlt\ 



q^\ + \ h{k + -) arctan( 



1 + it' 



Pk{q,t) 



,r2f e^' V^" ( ^0 



=7* 



/i \,(i + i^)' + ^o^' 



qm 



b. Hydrodynamic equations. The system of hydrodynamic equations for 
the velocity and density of the fluid 



dv dv 

-^ + iv + v—- = -e 

ot oq 

dp d , , ^ 



with specific initial conditions 



-It 



d_ 

dq 



27p dq' 



UJi 



+ -^e-'Y 



v{q, 0) = 0, pk{q, 0) = N^, exp ( -^g^ ) H^ 



n 



has solutions 

v{q,t) 

and pk{q,t) given by ([25]). 



-Wn 



1 + it, 



(l + |t)2+y,2^2^ 



c. r/ie potential Schrodinger-Burgers equation. The IVP for potential 
Schrodinger-Burgers equation 



dF 
'dt 



l,dF^ 



w, 



„2 



+ 7i^ + 7T(^^)' + Tr? 



ifi B'^F 



2'dq' 2^ 2 'dq^'' (26) 

^fc(g, 0) = t (f g2 _ ;, ln(iV,i7,(yf g))) , fc = 0, 1, 2, ... 



10 



has solutions Fk{q,t) = Fifc(g, t) + iF2,k{(l,t), where 



Fi, 



Un 



2 VI + i^. 



^TTWr^) '' + ^'"*^^ + '2^ -ctan(^), 



"2,fc 



Wo / i 

2 [{i + ity + wit^ 

-he-^^ In iffc 



g2 - he-^' In 



A^t 



-(7/2)t 



1/2- 



(1 + it)2 + wit^ 



iOo 



g7t 



2+2 



\ h \(i + ity + wit 



d. The Schrodinger-Burgers equation. The IVP for the non-hnear com- 
plex Schrodinger-Burgers equation 






^fc(g,o) = ? 



cJog — 2k\JhujQ 



Hk-iUfq) 



H, 



kU/fq) 



(27) 



has solutions 
Vk{q,t)-- 



-OOn 



Wo 



.1+i^ 



;i + it)2+«;2t^ 



;i + it)2 + wit 



27-2 



i2k\/huJoe 



-It 



3 7* 



1/2 



:i + ^t)2 + «;o2t2^ 



^fc-i(y^ ( (TTil^+^) ^) 



V 



H, 



k\\l h 1 nj_2i^24.,„2t2 ) 0') 



a I (1+24)2+^2*2 



.(28) 



The complex velocity function, written in the form V{q,t) = v{q,t) + 
iu{q,t), where v,u are real- valued functions, splits the Schrodinger-Burgers 
equation fl271) into the system 



dv du 



du dv dv h _^id'^v 

dt dq dq 2 dq'^ ' 



h ^..d'^u 
dq^ 



11 



This system with specific initial conditions 



HkLffq) 



v{q, 0) = 0, Uk{q, 0) = uqQ - 2kJhuJt 



clearly, has solutions v{q,t) and Uk{q,t) which are respectively, the real and 
imaginary parts of Vk{q,t) in expression ( !28|) . 

e. Motion of Zeros and Poles, Q^ = 0. From expression (125|) we see that 
the solution \l/fc(g, t) (also Pk{q,t) ) of the linear Schrodinger equation 
has zeros at points where 

H, 



q) = 0, 



\ h \{i + ^py + wit\ 

and these zeros are pole singularities of the solution Vk{q,t) (also |Vfc(g,t)p 
) for the nonlinear Schrodinger-Burgers equation (1271) . Therefore, for each 
A; = 0,l,2,3,..., the motion of the zeros and poles is described by the curves 

q^\t) = t!\1^ X e-^* X J {I + ^tf + wlt\ (29) 

where t^, \ I = l,2,...k, are the zeros of the Hermite polynomial Hk{^). 



Clearly, at initial time the position of the zeros and poles is g^ (0) = r^ Jh/uo, 

and when 7 > 0, t — > oo one has ql (t) — ?■ due to increasing mass /i(t) = e"^*, 
(dissipation). 

In the figure we illustrate the probability density function P2{q,t) for 
Caldirola-Kanai oscillator, which shows Dirac-delta behavior at time infinity 



and the position of the moving zeros and poles ^2 (^) = ""t? V ^"'^ ~'~ ^^^ ~'~ ^^' 

Q2 (^) = ~'l2 (^)- For the plots, constants are chosen as xq = h = uq = 1 
and 7 = 2. 
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Figure 1: Case Q^ = 0. a) Evolution of probability density P2{(l,t). b) Plot 
of moving zeros and poles. 

6 Conclusion 

In the present paper we have described time variable Madelung fluid and 
its linearization in terms of time variable linear Schrodinger equation. Our 
model, as descriptive of dissipative quantum fluid, admits exact solution for 
specific time dependent systems, like the harmonic oscillator with time de- 
pendent frequency and mass. In this case, exact time evolution has been 
described in terms of solutions for the corresponding damped classical oscil- 
lator. In particular, for the damping simulated by an exponentially grow- 
ing mass (the Caldirola-Kanai model), it can be shown that the quantum 
damping squeezes the density function of the fluid and leads to Dirac-delta 
function. This can have some implications in quantum cosmology. In fact, if 
\l/(g, t) is a solution of the Caldirola-Kanai oscillator for a damped system 



ih 



'dt 






-54-W-' 



then $(g,t) = '^*{q, — t), where (*) denotes complex conjugation, is a solu- 
tion of the related dual amplified system 






-7*^2 



g"*. 
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Hence, knowing that the solution \l/((3', t) of the Caldirola-Kanai model has 
merging zeros and describes collapse of the wave function to Dirac-delta func- 
tion at time infinity, leads to possible interpretation of the solution \l/*(g, —t) 
for the dual system. Namely, we will have expanding wave function with 
creation of zeros as point particles from initial singularity at time zero. And 
this evolution simulates quantum mechanism similar to creation of expanding 
Universe from initial singularity in Big-Bang cosmology. 

Finally, we note that by our results it is possible to find explicit exact 
solutions to a wide class of exactly solvable dissipative quantum fluids and 
complex Schrodinger-Burgers equations. For this we can use our recent work 
on exactly solvable dissipative systems, such as quantum Sturm-Liouville 
problems Ref.[7]. Then, it is possible also to describe the dynamics of the 
zeros and poles in the corresponding dissipative linear and nonlinear systems. 
These questions are under investigation now. 
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